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The anomalous Hall effect due to the surface conduction band of 3D topological insulators with an 
out-of-plane magnetization is always dominated by an intrinsic topological term of the order of the 
conductivity quantum. We determine the contributions due to the band structure, skew scattering, 
side jump and magnetic impurities on the same footing, demonstrating that the topological term, 
renormalized due to disorder, overwhelms all other terms, providing an unmistakable signature of 
Z2 topological order. Uncharacteristically, skew scattering contributes in the Born approximation 
as well as in the third order in the scattering potential, while in addition to the side-jump scattering 
term we identify a novel intrinsic side-jump term of a similar magnitude. These, however, never 
overwhelm the topological contribution. 



O 
> 

in 

m 



X 



I. INTRODUCTION 

Topological insulators (TIs) have been an intense focus 
of research and witnessed impressive experimental break- 
throughs in the past years 1 - — . Several materials were 
predicted to exhibit TI behavior—, and experiment has 
observed TI behavior in Bi Sb&, Bi 2 Se 3r ^ Bi 2 Te 3 ^^- 
and Sb2Teyli, with Heusler alloys also recently predicted 
to belong to this classic. This extraordinary novel state 
of matter is remarkable for being insulating in the bulk 
while conducting along the surface, as well as for the 
absence of backscattering, which leads to expected high 
mobilities, Klein tunneling, and the suppression of An- 
derson localization. 

The existence of chiral surface states is well estab- 
lished experimentally through angle-resolved photoe- 
mission (ARPES) and scanning-tunneling microscopy 
(STM) studies, yet in transport their contribution is 
overwhelmed by the bulk, which, due to unintentional 
defect-induced doping, is a low-density metal. A key 
aim of current physics research is the identification of 
a definitive signature of surface transport. One avenue is 
doping with magnetic impurities, whose moments align 
beyond a certain doping concentration and make TIs 
magnetic, opening a gap 1 -. Recent theories of mag- 
netic TIs include the Kerr and Faraday effects in a thin 
TI film coupled to a ferromagnetic and in a strong ex- 
ternal magnetic fielcU 5 -, spin transfer—, and anomalous 
magnetoresistanceii, and edge contributions 1 ^, while ex- 
periment has observed magnetization dynamics 19 and at- 
tempted to control surface magnetism 2 ^. The interplay 
of magnetism, strong spin-orbit coupling, disorder scat- 
tering and driving electric fields leads to the anomalous 
Hall effect (AHE) in magnetic TIs. The AHE has been 
a mainstay of condensed matter physics, studied exten- 
sively over the years^— , and the scattering mechanisms 
involved, skew scattering 2 ^, and side jum p 23 i 24 i 32 , are well 
known to be non-trivial. Furthermore, a fascinating in- 



stance of this effect has been predicted to occur in 2D 
magnetic TIs, where edge states are believed to give rise 
to a quantized AHE when the chemical potential lies 
in the magnetization-induced gap between the surface 
valence and conduction bands C 3 - (See Ref. [H for a de- 
scription of the transition between the ordinary and TI 
regimes in magnetotransport.) 

This work is devoted to the experimentally relevant 
problem of the AHE on the surface of a three-dimensional 
magnetic TI, focusing on the contribution from the sur- 
face conduction band. This already challenging prob- 
lem is complicated by the necessity of a quantum me- 
chanical starting point and of a matrix formulation de- 
scribing strong band structure spin-orbit coupling, spin- 
dependent scattering including relativistic corrections, 
Klein tunneling and Zitterbewegung on the same foot- 
ing. The fascinating and counterintuitive conclusion of 
our study is that a topological band structure contri- 
bution, known to exist in magnetic spin-orbit coupled 
system s 34 ' 35 , yields the dominant term, overwhelming all 
scattering contributions. The result, which is distinct 
from the topological magnetoelectric effect^ 3 -, is of the 
order of the conductivity quantum, independent of the 
magnetization, and unambiguously reflects Z2 topologi- 
cal order in transport. 

The outline of this paper is as follows. In Section II 
the Hamiltonian of the system is introduced. In Section 
III the kinetic equation for the density matrix is derived, 
containing all spin-dependent scattering terms in the first 
and second Born approximations. This includes the usual 
scalar scattering term in the first Born approximation, 
the skew scattering and side-jump scattering terms in 
the first Born approximation, and the usual skew scat- 
tering term in the second Born approximation (of third 
order in the scattering potential) . In Section IV the vari- 
ous nonequilibrium corrections to the density matrix are 
determined, and their contributions to the AHE are cal- 
culated. Comparison is made with existing experiments 
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as well as predictions for future experiments. Section V 
contains a summary and conclusions. 

II. HAMILTONIAN 

We study a topological insulator doped with mag- 
netic impurities, which give a net magnetic moment 
and spin-dependent scattering. The latter are described 
by H mag (r) = er ■ J2i^( r ~ Ri) s i, where the sum 
runs over the positions Ri of the magnetic ions with 
si, and the vector of Pauli spin matrices er represents 
the spins of the conduction electrons. We approximate 
V(r — Ri) = J S(r — Ri), with J the exchange con- 
stant between the localized moments and the itinerant 
carriers, and V the sample volume. The magnetic ions 
are assumed spin-polarized so that s/ = sz. Fourier 
transforming to the crystal momentum representation 
{|fe)}, with k the wave vector, the fe-diagonal term, 
•^mag = "mag J sa z = M <7 Z , gives the magnetization, 
with n mag the density of magnetic ions. The fc-off- 
diagonal term leads to spin-dependent scattering 

H *t'= ^£ ei(fe - fe,) - Rl - (1) 

We consider zero temperature, so that the splitting due 
to M is resolved, and the Fermi energy ep ^> |JVf| lies 
in the surface conduction band. We require Sft/H 3> 1 
for our theory to be applicable, with r the momentum 
relaxation time. The effective band Hamiltonian is 

H 0k = -Ak(T-8 + <T-M=^(T-in k , (2) 

where 9 is the tangential unit vector in 2D fc-space, and 
the constant A = 4.1eVA for f^Se^. Scalar terms are 
negligible at relevant transport densities 10 13 cm -2 — 
The eigenenergies are e± = ±V ' A 2 k 2 + M 2 . We will 
frequently use au = 2Ak/hQk and bk = 2M/f\£l k , so 
that a 2 + b\ = 1, and ap and bp at k = kp, with 
bp <C 1 at usual transport densities. Interaction with a 
static, uniform electric field E is contained in Hp kk i = 
^E kk' + Hpuui, the scalar part arising from the ordi- 
nary position operator Hp kk , = (eE ■ r) kk >l, with 1 
the identity matrix in spin space, and the side-jump part 
kk , = e\cr ■ (k x E) 5 kk / arising from the spin-orbit 

modification to the position operator, f — > f + Act x k 29 
with A a material-specific constant. Elastic scattering 
off charged impurities, static defects and magnetic ions 
(but not phonons or electrons) is contained in the scat- 
tering potential U kk > = U kk > J2j e i(fc ~ fe ' ) '' R, , with Rj 
the random locations of the impurities, which incorpo- 
rates the magnetic scattering term iJ^ag °f ^1- ©■ 
The potential of a single magnetic impurity U kk > — 
(1 — iXcr ■ k x k')U kk > — (J s /V) a z , with U kk i the ma- 
trix element between plane waves, and Xkp 1, while 
J = for non-magnetic impurities. The full Hamiltonian 



H k = H ok + Hpkk' + U kk < + -ff^ag - We do not include 
weak localization, which in the regime epr/h 3> 1 pro- 
duce only a small correction to the conductivity, and is 
known to be unimportant in the AHEi^ 

It is well known that in magnetic TIs the Berry cur- 
vatures of the surface valence and conduction bands are 
different, so that when the chemical potential is in the 
band gap the conductivity is e 2 /2/t. 35 ' 37 This quantity is 
regarded as a constant offset in our work, which is con- 
cerned rather with the response of the conduction elec- 
trons, that is, the surface conduction band. 

Nontrivial processes such as spin-flip scattering and the 
Kondo effect can be accounted for in our theory, provided 
that one works in the grand canonical ensemble in order 
to account fully for the spin-flip factors. Several reasons 
compel us to avoid doing so. Given that the magnetiza- 
tion must be large enough to be resolved, the local mag- 
netic moments are assumed to be all lined up. Moreover, 
although it is reasonable to assume that the Kondo effect 
may affect the diagonal term in the conductivity signif- 
icantly, there is no evidence that it should influence the 
AHE. We will find in this work that all corrections to the 
AHE are zeroth order or higher in the scattering poten- 
tial, in other words next-to-leading order. Considering 
the very few studies of the Kondo effect in topological 
insulators^ we have no clear indication of how high the 
Kondo temperature is expected to be, and the nontrivial 
effect of spin-momentum locking on Kondo singlet for- 
mation has not been determined. These factors will be 
considered in a forthcoming paper. 



III. KINETIC EQUATION 

The Liouville equation for the density operator p is 
projected onto the basis {|fe)} as in Ref. |36|. The density 
matrix p kk > = f k <5 fcfc / + g kk > , where g kk > is off-diagonal 
in fc, and all quantities are matrices in spin space. To 
first order in E, the diagonal part f k satisfies 3 ^ 

< ^+ l j i [Hv k J k ] + J{f k )=V k . (3) 

The scattering term J(f k ) is given below, the driv- 
ing term V k = -^[H^,f ok \, and f ok is the equilib- 
rium density matrix, which is diagonal in k. We write 
fk = ri fcl + Sk, where Sk is a 2 x 2 Hcrmitian matrix. 
Every matrix in this problem can be written in terms of 
a scalar part, labeled by the subscript n, and er. Rather 
than choosing Cartesian coordinates to express the latter, 
it is more natural to identify three orthogonal directions 
in reciprocal space, which we denote by tt k , k, and z e ff, 
and project cr onto them. The three orthogonal direc- 
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tions in reciprocal space that we use for reference are 
flk = — ak + b k z 
k eff = k 



(4) 



z e ff = a k z + b k 6. 



We project er as <J k ,\\ = n ' ^k, &k,k = cr ■ k, and 
o'fe.ze/f = ' ^eff- Note that <7fe n commutes with -f/ofc, 
while (Jk,k and <J k>Ze tj do not. Therefore we shall often 
use the abbreviation _L to refer to vectors in the plane 
spanned by k and z e ff. We project S k onto the three 



directions in fc-space, writing Sk — S k n 
and defining S ht \\ = (1/2) s ki \\ a k ^\, S k , k 
and Sk, Zeff = (1/ 2) Sfe,* e// 0"fe,*e// ■ 
Equation (|18ap can be written as 



: (1/2) S k;k &k,k 



dn k 
dS k \\ 



dS k ._ 
dt 



PnJ(fk) = X>fe,n, (5a) 
P\\J(fk) = ©fej, (5b) 
[H k ,S kA ]+Pj{f k ) = D fe)± . (5c) 



The projector P|| acts on a matrix .M as tr (M.<J k .\\ ), 
with tr the matrix trace, while Pj_ singles out the part 
orthogonal to Ho k . The total scattering term J(f k ) = 
J Born (fk) + J 3rd (fk), with 



J 



Vk) 



dt' ~ 



kk - 



(6) 



Considering the negligible size of b k at k = kp, the con- 
dition yielding the suppression of backscattering is effec- 
tively unmodified from the non-magnetic case. In the 
absence of scalar terms in the Hamiltonian the scalar 
and spin-dependent parts of the density matrix, n k and 
S k , are decoupled in the first Born approximation. In 
the second Born approximation we obtain the additional 
scattering term 



J 6r Vk) 



dt' 



dt"[U, 



[U,e 



-[uj}< 



kk, 



(7) 



while () represents averaging over impurity configura- 
tions. The former can be further broken down into 
J Born (h) = Jo(fk) + Jfs 0rn (fk) + J Born (fk). The first 
term, Jo(f k ), in which A = E = in the time evolution 
operator, represents elastic, spin-independent, pure mo- 
mentum scattering. In Jf° rn {f k ) we allow A to be finite 
but E = 0. In Jf° rn {f k ) both A and E are finite, thus 

J Born (fk) acts on the equilibrium density matrix fo k . 
Beyond the Born approximation we retain the leading 
term J 3rd {fa) = J^ d (f k ), with A finite but E = 0, and 
which is customarily responsible for skew scattering™ 
The contribution due to magnetic impurities is also oc er 2 



and is contained in jf s orn (f k ). 

The formalism we adopt has a quantum mechanical 
foundation and accounts for the three mechanisms known 
to be important in the AHE. The interplay of band struc- 
ture spin-orbit interaction, adiabatic change in k in the 
external electric field, and out-of-plane magnetization 
leads to a sideways displacement of carriers even in the 
absence of scattering {intrinsic)^. The spin dependence 
of the impurity potentials via exchange (J) and extrin- 
sic spin-orbit coupling causes asymmetric scattering of 
up and down spins (skew scattering)^, and a sideways 
displacement during scattering (side jump)22. 



A. Spin-independent scattering term in the Born approximation 

The projections of Jo(f k ) that we will need in this work, derived in an analogous manner to Ref. l36l . are 
Tin- f d 2 k' 

p \\MSk\\) = j j^y\U kk '\ 2 {s kl \~s k , ll ){l + bl+a 2 k cos-f)a kll S{e' + -e + ) 

" 7TTI ' f fc' " A A 

P±Jo(S k \\) = — jjfi" J (2^)2 \ u kk'\ 2 (s k \\ - Sfe'ii) [(n fc / ■ k)o k , k + (fi fe / • z e ff)a k ,z e ff] S(e' + - e+) (8) 
nn f d 2 k' 

P\\Jo(S k ±) =--^r / \U kk ,\ 2 [(s k , eff + s' keff )(Cl k , ■ k) + (s Zieff - s' zeff ) (fi fe , ■ z eff )] a k]l 8(s' + - s+), 



4 



where rii is the impurity density, which emerges once one averages over impurity configurations. 



B. Born approximation skew scattering and magnetic impurity scattering 

The effect of this term on the scalar part of the density matrix gives a scalar contribution, which does not yield an 
electrical current. Next, Jf s orn acting on the spin-dependent part of the density matrix gives 

J* orn {S k ) = J |^ \U kk> \ 2 sin 7 [ ( T ze -^' t '(5 fe - S w )e iH ** e"^"' (a z S k S kl a z )e iH ^'] + h.c. (9) 

In this subsection we omit factors of ft in the time evolution operators, which are restored in the final results. The 
two terms in the square brackets are divided as J^ orn (S k ) = Jg°[ n (S k ) + Jg ° 2 r "(5 fc ). The first yields 

7TT? \k 2 C d 2 k' 

JssT^) = —5tT (T '] J2vf I^'| 2 sin7(5 fc - S k ,) ■ (fi fe x fi fe ,) S(e> + e+ ). (10) 

This scattering term becomes a source term in the kinetic equation. One will only obtain a contribution linear in r 
if one feeds S k \\ into the RHS, since S k \\ oc t. But at the same time S k \\ cx fife, and its contribution vanishes. The 
surviving scattering term is j^°™(Sfc), and we refer to this term simply as Jf s orn (S k ). Its parallel projection is 



p \\J?s° rn ( S k) = mi ^ h bk a k,\\ J |Wfcfc'| 2 sin7a fc [sin 7 (s fez , e// - s k ' z ,eff) - M 1 ~ cosj) (s k - s' k )] S(e' + - e+) 

imi\k 2 a k f d 2 k' 2 . 

2^ Cfe,|| J j2^\ Ukk '\ Sln 7{(Sfe,e// + S '=',e//)(l + COS7) 

- Sin7 [a fc (.Sfe|| - s k >\\) - h(skz,eff - Sfe'z.e//)]} S(e' + - £+). 

(11) 

It does not contribute to the Hall current. Its perpendicular projections are 
TTn,iXk 2 ak f d 2 k' 

PkJf s orn {S k ) = — % — crfc, fe J 7 ^ 2 -|Wfefe/| 2 sin7[a fe sin7(s fe + 4)- (1 - COS7) {s k \\ - s k >\\)\ S(e' + - e+) 

irn\k 2 a 2 f d 2 kl 

P z ,effJf s ° rn {S k ) = k a k , Zeff J -— ^ \U kk ,\ 2 sin 7 [2 sin 7 (s z , eff - s z , eff ,) - b k (l - cos 7 ) (s k - s' k )\ 5(e' + - e+) 



+ 



(27T) 

77ll,\k 2 (li, '' 1,2 



/d k 
—-2 \U k k> | 2 sin 2 7 b k (sfe|| - Sfc,|| ) <5(e' + - £+). 



2ft 

.... (12) 

Next, in exactly the same manner as above, taking into account the magnetic term in U kk i cx J, we obtain the 
scattering term due to magnetic impurities to linear order in bp (i.e. linear in J), referred to as Jma g n (fk)- We find 

that Jmaq n {Sk) is a scalar and therefore does not contribute to transport, while 



iBorm \ n mag Jsa k kU eErcosO df 0+ 

'mag \ n k) = ~ ' 

where U = {l/2w) J djU kk '('y) and n mag < m. 



jBomt \ n mag jsa k Ku em cos u a j 0+ „ n „, 
J mag (n*) = ^ 2ft dk a ' Zen ' { 3) 



C. Skew scattering beyond the Born approximation 



The skew scattering term of third order in the scattering potential is 

Jll d (h) = </ dt ' e ~"' J dt " ^ ft, ^ lHat [U, e- lHat [U, f] e lHat ] e^ 4 ]„„). (14) 

In each term the A-dependent part can appear in three places, so we have a total of 12 terms to evaluate. In each term 
in we retain the imaginary, spin-dependent part. We find that only J^l d (n k ) contributes to skew scattering, so 
we only retain this term to simplify the algebra below. We let k\ , ki — > k 2 since this term will eventually be evaluated 
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at the Fermi energy. Skipping a large amount of laborious, yet straightforward, detail, the parallel projection of the 
third order skew scattering term is 



Xk 2 b k 

Pl\J*I d = —nr- \U\ 3 (n k2 ~ n kl ) (cr ■ fi fe ) [ (9 + 3b 2 k - a 2 k )(smj kkl + smj klk2 + sm^k) 



4ft 



3a 2 

(sin27 fcfcl + sin27 fclfc2 + sin27 fc2fc )] D klk D k2 k- 



(15) 



This term vanishes when we integrate over k\ and k 2 , thus there is no skew scattering contribution to the AHE that 
is oc r. The projections || k can be written as 



P k J: 



3rd. A 



2a k b k Xk 2 
ft 3 

2a k b k \k 2 
ft 3 





( n k 2 


- n k ) (cr • 


fc)[- 


- sin 2 7 fcfcl 


-1- sin7fc lfc2 sin7 fcfcl 


- sin7 fcfcl sin7 fefc2 ] 


M 3 


(n kl 


- rife) (cr • 


fc)[- 


- sin 2 7fe fe2 


- sin7fc lfc2 sin7 fcfc2 


- sin7 fcfcl sin7 fcfc2 ] 


M 3 


( n k 2 


- rifcj (cr 


■fc)( 


- sin 2 7 fcfcl 


+ sin7 fclfc2 sin7 fcfcl 


- sin7fefe 1 sin7 fcfc2 ) 


M 3 


{n k2 


- nfej (<T 


■fc)( 


sin 2 7fcfc 2 + 


sin7 fcl fc 2 sin7fefe 2 + 


sin7fc fcl sin7fcfc 2 ). 



(16) 



h3 

1aub k \k 2 
+ P 

PkJss vanishes for isotropic scattering and does not contribute to the Hall current. The projections || z e ff can be 
written as 



Pz,effJ. 



P z,effJ, 



Xk^d 

! - \U\ 3 {n k2 - rife) (cr • £eff) { sin7fefe 1 [(l + 6 2 )(1 + cos7 fcfcl ) + a 2 k cosj kk2 + a 2 k cos^ klk2 ] 

+ sin7fe lfe2 [(l + 6 2 )(1 - cos7fc fcl ) + a 2 k cosj kk2 - a 2 k cosj klk2 ] 
+ sin7fe 2 fe[(l + 6 2 )(1 - cos7 fcfcl ) - a 2 k cosj kk . 2 + a 2 k cos-f klk2 }} + (fei O fc 2 ) 

l^l 3 ( n k 2 - nki) ( CT ' *e//) { shi7fefei[(l + & fe)(l + cos7fc fcl ) + a 2 k cos~/ kk . 2 + a 2 k cosj klk2 ] 

+ shi7fe lfe2 [(l + 6 2 )(1 - cos7fc fcl ) + a 2 cos7fe fe2 - a 2 k cosj klk2 ] 
+ sir±7fe 2 fe[(l + 6 2 )(1 - cos7 fcfcl ) - a 2 cos 7^ + a\ cos7 fclfc ,]} + {k\ o k 2 ) 

, 01) 

I 



3rd,B 



Xk 2 a k 



We will leave the results in this form until Section IV, In the Born approximation, we solve Eq. (|18b[) to first 
where they will be evaluated explicitly. order in U 

D. Side jump 



One contribution due to side jump arises from the in- 
trinsic side-jump driving term, considered in Section IV. 
The other two contributions come from corrections to 
the scattering term. We retain only the parallel projec- 
tion of the side-jump scattering terms, the others being 
of higher order in H/(bft). The easiest way is to write 
the kinetic equation in terms of the k— diagonal and off- 
diagonal components of the density matrix 

^ + ~ [H ok J k ] = - % - [H*;, / fe ] - I [UMWa) 



(18b) 



dt' e 



-iHf 



U,f(t-t') 



JHt' 



(19) 



fefe' 



To include side jump, we take into account an additional 
correction to g arising from He- Still in the Born ap- 
proximation we write g = go + gE 



6 



dt h l U '-"J"- ^ 



| ^ dt 7 e-"*' e- ii?0 *' [17, /(t - t')] e^ *' 



dg E ,kk> , i ,{r - 1 * rrr ~ i 
(ft ft [-"QjfiBjfcfe' — ~T 9oJfcfc' 



fcfc ' (20) 



9EMk' = - ~ e- i6ot '[H E ,go(t - t')} e^ ot '| 



kk' 



Factors of h have been left out of the time evolution operators to shorten the equations. We adopt the Markovian 
approximation f(t — t') ~ f(t) = /, and, remembering that He is first order in the electric field, we replace / — > /o, 
the equilibrium density operator 

9E.kk>{t)~-^{ r dfe-^'U-^^HE^-^^Ujo]^^']^^''] . (21) 
6 Uo J fcfc' 

Both /fg and t/ have scalar and spin-orbit parts. We first evaluate the easier contribution from H^ and the scalar 
part of the scattering potential, which we call g^ k k ■ This contribution can be expressed as 

9 S lk lk2 = — h I dt' e-* e- iH ^ (H% ki g 0klk2 - Qok^J e iH ^ . (22) 
where gok x k 2 is taken to zeroth order in A. This gives rise to the scattering term 

Jsj(fok) = \Y1 Uk kl g s i Mk + h - c > 



(23) 



where only the scalar part of Ukk 1 contributes, since g^ kik is already first order in A. 
We require in addition the scattering term 



Jsc(fok) = r U ^x9E,kxk + h - c - (24) 



Ti- 
ki 



The term J sc {fok) has two contributions. One contribution arises from the spin-independent part of |J7fefe'| and go to 
first order in A. The other term comes from the spin-dependent part of |f/fefe/| and go to zeroth order in A. Therefore 
we only need to evaluate one other term to first order in A, 

if 00 ,, 

9E M k 3 =-JiJ o dt ' e ^ e ~ lH ° hlt ( H Ek lk2 90k 2 k 3 - gok lk2 H E c k2k3 ) e^-s* . (25) 
We will leave the result in this form until Section IV when i t will be evaluated explicitly. 



IV. ANOMALOUS HALL EFFECT 

We solve Eq. ([5]) perturbatively in the small parameter 
}i/(sft) (which is proportional to the impurity density 
rii) as described in Ref. .36, as well as in A. In transport 
this expansion starts at order [h/ (eft)] -1 (i.e. ), 
reflecting the competition between the driving electric 
field and impurity scattering resulting in a shift of the 
Fermi surface. The leading terms in [h/ '(e^r)] -1 appear 
in rik and Sfc,||— , namely in the components of the density 
matrix which commute with Hok- Identifying these terms 



correctly is crucial, since if terms cx n i exist in the AHE 
they will be dominant at high mobilities. 

We first solve Eqs. ([5]) to zeroth order in A, setting 
J(fk) — > Jo(fk)- We refer to this part of the solution 
as the bare part, which is gives the intrinsic contribu- 
tion to the AHE (i.e. A = 0.) Once the bare compo- 
nents of Sk are determined, we project repeatedly be- 
tween Sfeii and Sk,± in Eq. ([5]) until we have found all 
terms to ^/(e^r)] -1 and [ft/(e f r)]°. These projections 
give disorder renormalizations, explained in a separate 
subsection below. Subsequently, the scattering terms 
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Jf s orn (fk) + J 3rd {fk) act on the bare solution to gen- 
erate new driving terms to first order in A, for which an 
analogous solution is found to the density matrix. These 
give the skew scattering contributions, and the contribu- 
tion due to magnetic impurities. Finally the side-jump 
driving terms are found, in which f k — > Jofc- Given the 
lengthy expressions, below we present only the terms that 
contribute to the AHE. 

The appearance of terms to order zero in K/(eft) is 
standard in transport and important in spin-dependent 
transport. These terms depend on the angular structure 
of the scattering potential, thus it is not enough to study 
only the band structure contribution, but also its disorder 
renormalization. Disorder renormalizations are crucial 
in spin-related transport, where they can cancel band 
structure contributions^. Though additional terms of 
order [H/(eFT))° may exist, they arise from scattering 
terms oc nf , describing quantum interference and related 
effects relevant at higher impurity densities. 



A. Bare driving term 

To first order in E, we write fk = fok + fEk, with cor- 
responding notation for nk and Sk- The Fermi-Dirac 
function for the positive energy branch is denoted by 
/o+ = /o(e+), with /o_ = for electron doping. The 
equilibrium density matrix is fo k = \ (/o+ + /o-)l + 
\ (/o+ - /o— )crfe,||- The driving term is 



The transport scattering time for screened Coulomb im- 
purities is defined as 



1 _ £+ n t W k 
t~ A 2HA 

Z 2 e i 



(Co ~ Ci) 



(29) 



bt 



a\ COS7)- 



k T F ' 

2k F / 



C contains a factor of 1/2 as it is defined in analogy with 
the B = case^ while (0,1 refer to its Fourier compo- 
nents, and = ?f, with the Wigner-Seitz radius r s — 

e 2 /{2^e t r A). This way \U ky \ 2 = W k /(wa.% + r -ff. 

The (bare) perpendicular part of the spin density ma- 
trix is 



>Ek.. 



1 

2° 



Z eff 



eE-6 

hkfl 



o-k fo-\ 



(30) 



0+- 



This completes the evaluation of the bare correction to 
the spin density matrix, n b £ k e + S^ k e , which will yield the 
intrinsic contribution to the Hall current and its disorder 
renormalization. 



Disorder renormalizations 



Vi. 



[H k , fok] 



= {eE / H ).^t- i -^[a.(kxE)J 0k ]. 

Its projections are V k = V k . n + V kt \\ + V k M + T^k 
To zeroth order in A 



(26) 



V, 



eE-k f df Q+ 

1 eE-k/df 0+ 
2 ah ^~^{~dk~ 

eE6 



(27) 



1 eE6 

T^k,k = - &k,k — 7~j (Ik (J0+) 

2 hk 



_ 1 eE-k 

U KZeJf - ~ O U k,Z ef } fcJ , a kO k (JO+J- 



hk 



Henceforth E II x. We obtain 



bare 
l Ek 



eEr ■ k d/ 0H 



2h 



dk 



Ek,\\ 



1 - eEr ■ k df 0+ 

-*n k — — — 



(28) 



The correction to the bare density matrix is renormal- 
ized by scattering between S k ^± and S k ti. In addition, a 
number of terms appearing below also need to be renor- 
malized in this way. We briefly describe the technical 
procedure in this subsection. 

Projections of the form — P± J(S kz \\) involve no compli- 
cations, since Eq. (|5c[) is solved immediately by applying 
the time evolution operator. 3 - For projections of the form 
—P\\J(Sk.±), it is helpful to observe that the Fourier ex- 
pansions of Sk and s z . e ff throughout this work will only 
have Fourier component 1 (i.e. e ±ie ), thus generally, us- 
ing Eqs., 



u t / q \ a kbkSk 

-P\\Jo(JEk,k) = — Ofcji 



" \ Q-kbkSzeff 
-P\\Jo{SEk,z cl} ) — — cr fej 



(31) 



The effective scattering times r M , r c+ are given below. 

The projected terms —P±Jo(S k ,\\) an d —P\\J(S k ,±) then 
act as new driving terms for Sk,± and Sk,\\ respectively, in 
effect renormalizing the bare results by a multiplicative 
constant which depends on the angular characteristics 
of the scattering potential. All terms in this problem 
have to be renormalized in this way. Yet all terms in 
the density matrix consist of Fourier component 1 (that 
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is, they are of the form e ±n , or alternatively sin 7 and 
cos 7), thus the renormalizations all follow the pattern of 
Eqs. ((51)1 . 

Only the bare term S^ e needs to be renormalizcd 
due to scattering from 5 fe || to Sk±. The corrections to 
the perpendicular part of the spin density matrix due to 
scattering from Sg£u into Sk,k an d Sk, Zcff are 



QVtX 

D Ek,k 



eE x b k r cos df 



?iT c _f2 dk 
eE x r sin df 0+ 
fi.T s f2 dk 

The effective scattering times are 



• k 



(32) 



S 



Ek.. 



Z eff, 



1 
1 
1 

0(7) 
s(7) 

m(t) 



WT (co±Cl) 

1 cos 7 

J(sinf + f) 2 
1 sm7 

^ (sini + ^) 2 
1 1 — cos 7 

2( S ini + ^) 2 ' 



(33) 



The renormalizations of all the other terms appearing be- 
low due to scattering from S k \\ to Sk± are of higher orders 
in A or hj {epr) and need not be considered. Renormal- 
izations due to scattering from Sk± to Sk\\ can be worked 
out immediately from Eq. (|8]). 

We emphasize that, aside from the transport relaxation 
time t, the effective scattering times appearing in this 
work are defined for convenience and they do not nec- 
essarily represent specific scattering processes, but are 
auxiliary quantities in achieving a final result. 



C. Extrinsic contributions to the density matrix 



full spin-dependent Ukk' in both, generating new source 
terms oc A in the kinetic equation. The terms linear in A 
in J Born (f E k) yield 



q 



Ek,k 



2 eE x cos9 t - 
-Afc a k b k — w — d(k- k F )cr ■ k. 



Ml 



(34) 



The effective scattering time Tf s orn is defined as 



2HA 2 



(vo - m) 



sin 2 7 



(35) 



( sin ? 



Magnetic impurity scattering is also part of J rn (fEk)- 
The scattering term Jma g n ( n k) evaluated in Eq. (fi"3)) 
yields a driving term for SEk±i which gives a correction 
cx n,;T and therefore independent of the impurity density 



qmag 

°Ek± — "fc 



r eE cos 9 



1 mag 



S(k — kp) <x ■ k. 



(36) 



Here l/r mag = n mag JskU/(2h 2 A). 

We come to the skew scattering term beyond the Born 
approximation, JsI d (fEk)- In the term J^l d {f E a k e ) om y 
JsI d ( nb Ek e ) i s rea l an d finite, and we find that the pro- 
jection PkJgI d does not contribute to the AHE. We 
concentrate on the sum of the terms P z , e ffJsI d ' A + 

Pz,effJ!s d ' B = p z,effJfs d - This scattering term will act 
on the electric-field induced correction to the scalar part 
of the density matrix, n E k, which has the angular struc- 
ture riEk = n{k) cos9. In light of this detail, we can boil 
the term P z , e ffJ^l d to 



P^effJlf = \k 2 a k n(k) r^l (a ■ z eff ) 



suit 

ss 



(37) 



r 3rd 



2hA A 



The next step comprises two parts. Firstly, we feed 
into J Born {fk) and J 3rd (/ fc ), with the 



-.bare 
L Ek 



Qbare 
°Ek 



where the following auxiliary quantities have been de- 
fined 



1 f f d9 l f d9 2 F s z sin 0(2 cos 2 -cos0-cos0i) 

nj J 2^ J 2^ ( si n^ + ^)(sin^ + ^)(sinlfi + ^) 

sin7 fcfcl [(l + cos7 fcfcl ) + cos7 fcfc2 + cos7 fclfc2 ] + shi7 fclfe2 [(l - cos7 fcfcl ) + cos7 fcfc2 - cos7 fclfc2 ] 
+ sin7fc 2 fc[(l - coS7fefe 1 ) - cos7 fcfc2 + coS7fe 1 fe 2 ]. 

I 




The contribution of this scattering source term to the spin density matrix is 

, S s,3rd _ Xk 2 ak_n(k) sin 



r,ss,3rd ™ °' 1LU / £,\ fon\ 
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The term S k ' k 3rd found in Eq. (|39)) above does not con- 
tribute to the Hall current, but its disorder renormaliza- 
tion yields, 



riss,3rd 11 2 2 7 



k°k 



eE sin( 



■ 5(k - hp) (T fc> ||. 



(40) 

This term contributes to the Hall current. It is the only 
skew scattering contribution to the anomalous Hall effect, 
and its value is determined explicitly below. 



Finally, the side-jump scattering term J^°™(/ofc) con- 
tains both electric field interaction terms H E c kk , +He kk , 
in the time evolution operators. It is oc [h/ '(e^r)] -1 , 
thus only its parallel projection yields a term of order 
[h/ (epr)] . The side-jump scattering term consists of 
J"e (fok) + Jsc(fok), defined in Eqs. j23J and (gU respec- 
tively. Substituting for npk we obtain for Pu Jg (/ofe), 
after a lot of algebra, 



P\\J S E j (f0k) 



nnieb k \\E\ df + 



d 2 k' 



\U kk ,\ 2 (1 + fl k ■ flu) (k y - k') 5(e' - e+). 



h de+ J (2tt) 

The term J sc (fok) has two contributions: one from the spin-independent part of |f/fefe'| and go to first order in A 



(41) 



P\\Jsc(f0k) 



nebX\E\ df ok f d 2 k' 
■ &k, 



h ~ K '" de k J (2tt) 2 
and another from the spin-dependent part of |t/fefe'| and go to zeroth order in A 



\U kk ' | 2 {k y - k' y ) 5(s k < - £fc), 



P^Ufok) 



nriieXbk 2 dfo+ 
2h ~dsZ 



d 2 k' - 

/r( x ? \U kk > \ 2 sin 7 fcfc / (1 + fl k , ■ flfc) 



<5(e+ - fn 



(42) 



(43) 



r 



Integrating 

P\\JUh 



by parts and adding P\\J} c (fok) + 
o k ) — P\\JE(fok) yields an overall factor of 
2 in the parallel projection of the side-jump scattering 
term. The overall result for the side-jump scattering 
contribution is 



qJsj _ 
Ek,\\ - 



Aeb k XE x e + 
A 2 



sin 6 5{k — kp) o k 



(44) 



A further contribution due to side jump arises from the 
intrinsic side-jump driving term, 



T [H E , fok] 



ea k EXk v 



(f 0+ -fo-)<T-k. (45) 



This term only contributes to the perpendicular driving 
term and yields a correction to SE k ± , found straightfor- 
wardly evaluated by applying the time evolution opera- 
tor. We note in passing that this term is also responsible 
for the survival of side-jump spin-Hall effect in systems 
with band structure spin-orbit coupling 4 ^-. In the AHE 
in TI it yields 



S 



sj,int 
Ek,z c f 



ea k E x Xk sin 9 

m 



fo+ 



(46) 



This term is projected back onto S k \\ as described in Sec. 
HVB| giving an additional renormalization. Both Eqs. 
(|44| and (|46|) contribute to the AHE. 

The current operator j has contributions from the 



band Hamiltonian, jo 



er x z, as well as j E 



2£\\ 

h 



<x x E, and j v = er x (fc - k')U k k>- These 
latter two cancel, as they represent the force acting on 
the system. The central result of our work is the AHE 
conductivity, which, ignoring terms of order b 2 F , can be 
divided into two terms (int = intrinsic, ext = extrinsic) 



e 

2h 



2h 



b F (Xkjr) 9 



bpT 



St 



(t c _ + T s 


k t ] 




/ ' mag . 


T 


T2 ) 


^ T Born ' 
ss 





(47) 

All rs are defined in the supplement, and 1/T mag on bp, 
rendering this term negligible. Both Sf£™ and S k a k e 

contribute — |^ to oj™ . The total (— can also be ex- 
pressed in terms of the Berry curvature^, and is thus 
a topological quantity. It is similar to the result of 
Ref. |35j in the vicinity of a ferromagnetic layer, and can 
be identified with a monopole. The corrections are dis- 
order renormalizations due to scattering between S k ^± 
and S kt \\. In cr^f the first two terms in brackets rep- 
resent the combined contribution of side jump scatter- 
ing and intrinsic side jump plus their disorder renormal- 
izations, and the last two terms represent the contribu- 
tions due to skew scattering in the Born approximation 
and beyond it respectively. An important unprecedented 
finding of our study is that neither extrinsic spin-orbit 
coupling nor magnetic impurity scattering give a driv- 
ing term in the kinetic equation parallel to Ho k (that is, 
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|| Jlfe), thus the AHE response does not contain a term 
of order [fi,/(£^r)] _1 due to extrinsic spin-orbit scatter- 
ing or magnetic impurity scattering. Such a term would 
have overwhelmed the intrinsic topological term in the 
ballistic limit. Since the intrinsic topological term is also 
of order [H/(£ft)]° , we conclude that there is no term 
oc n" 1 (i.e. oc r only) in the AHE response of TI. 

In TI the spin and charge degrees of freedom are inher- 
ently coupled and the AHE current can also be viewed 
a steady-state in-plane spin polarization in a direction 
parallel to the electric field. The (bare) topological term 
has two equal contributions, which are part of the cor- 
rection to the density matrix orthogonal to the effective 
Zeeman field, therefore they represent an electric-field in- 
duced displacement of the spin in a direction transverse 
to its original direction. They are also obtained in the 
Heisenberg equation of motion if one takes into account 
the fact that k is changing adiabatically. Physically, the 
cc-component of the effective Zeeman field is chang- 
ing adiabatically, and the out-of-plane spin component 
undergoes a small rotation about this new effective field. 
Consequently, each spin acquires a steady state compo- 
nent parallel to E, which in turn causes k to acquire 
a small component in the direction perpendicular to E. 
Elastic, pure momentum scattering [contained in Jo(fk)] 
reduces this spin polarization because, in scattering from 
one point on the Fermi surface to another, the spin has 
to line up with a different effective field f2fc. The extra 
spin component of each electron is oc M, however the 
final result is independent of M, as the integrand con- 
tains a monopole located at the origin in fc-spac o 34 ' 35 . As 
M the effect disappears, since the correction to the 
orthogonal part of the density matrix vanishes. Finally, 
though the form of this term would hint that it is observ- 
able for infinitcsimally small M, it was tacitly assumed 
that M exceeds the disorder and thermal broadening. 

Experimentally, for charged impurity scattering the 
figures depend on the Wigner-Seitz radius r s , represent- 
ing the ratio of the Coulomb interaction energy and the 
kinetic energy, with e r the relative permittivity. For 
Bi 2 Se 3 with r s — 0.142&, the dominant term by far is 



a l y n J w -0.53 (e 2 /2h) w -e 2 /4h. We can consider also 
the (artificial) limit r s — > 0, which implies e r — > oo, arti- 
ficially turning off the Coulomb interaction. As r s — > 0, 
the prefactor of — e 2 /2h tends to 0.61, while at r s = 4, 
the limit of RPA in this case, it is s» 0.12. Interestingly, 
for short-range scattering, the AHE current changes sign, 
with cr™* w 0.18 (e 2 /2h). We note that A is not known 
for TI, yet that does not affect the central result of this 
work, since the intrinsic term is dominant. We have not 
taken into account the term in Jfj d linear in E, compa- 
rable in magnitude to the side jump. 

We recall that the result presented above represents 
the contribution from the conduction band, and in prin- 
ciple an extra e 2 /2h needs to be added to the total result 
to obtain the signal expected in experiment. Since the 
conduction band is expected to contribute s» — e 2 /4/i, 
this does not make a difference in absolute terms. Com- 
parison with experiment must await further advances in 
materials growth, yet the most recent developments offer 
renewed grounds for hope.— 



V. SUMMARY 

We have determined the AHE conductivity due to the 
surface conduction band of 3D TI, including all intrin- 
sic and extrinsic contributions. We have identified a 
skew scattering term in the Born approximation and 
one of third order in the scattering potential, as well 
as an intrinsic side-jump term and a side jump scat- 
tering term. The intrinsic topological term, renormal- 
ized by disorder, is dominant and of the order of e 2 /4/i, 
though non-universal, and easily observable experimen- 
tally. This finding provides an unmistakable signature of 
surface transport in TI. 
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